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ABSTRACT. Asymptotic decay laws for planar and nonplanar shock waves and the first 
order associated discontinuities that catch up with the shock from behind are obtained using 
four different approximation methods. The singular surface theory is used to derive a pair 
of transport equations for the shock strength and the associated first order discontinuity, 
which represents the effect of precursor disturbances that overtake the shock from behind. 
The asymptotic behaviour of both the discontinuities is completely analysed. It is noticed 
that the decay of a first order discontinuity is much faster than the decay of the shock; 
indeed, if the amplitude of the accompanying discontinuity is small then the shock decays 
faster as compared to the case when the amplitude of the first order discontinuity is finite 
(not necessarily small). It is shown that for a weak shock, the precursor disturbance evolves 
like an acceleration wave at the leading order. We show that the asymptotic decay laws for 
weak shocks and the accompanying first order discontinuity are exactly the ones obtained 
by using the theory of nonlinear geometrical optics, the theory of simple waves using 
Riemann invariants, and the theory of relatively undistorted waves. It follows that the 
relatively undistorted wave approximation is a consequence of the simple wave formalism 
using Riemann invariants. 



1. Introduction 

The study of the evolutionary behaviour of nonlinear waves in diverse branches of con- 
tinuum mechanics has long been a subject of great interest from both mathematical and 
physical points of view. A rigorous mathematical approach to describe the kinematics of a 
shock has been proposed by Maslov [ 1 1 using the theory of generalized functions. General 
properties and the instantaneous growth and decay behaviour of shock waves in elastic 
materials have been studied using the singular surface theory in a series of articles by Chen 
[2|, Nunziato and Herrmann , Chen and Gurtin [4], Eringen and Suhubi |5|, McCarthy 
[6 1, Ting [7 1, Bailey and Chen [8| and, Quintanilla and Straughan [9]. Evolutionary be- 
haviour of shocks in fluids, using singular surface theory, has been discussed by Grinfeld 
iflOl Anile El, Straughan fl2] Jordan ifMllBl . Radha et alEQ, and Pandey et al 
fPTl . Using a procedure based on the kinematics of one-dimensional motion, Sharma and 
Radha lfl8l [191 and Batt and Ravindran ||20l have studied the evolutionary behaviour of 
shocks. Numerical approaches such as the Godunov method have also been used to under- 
stand shock evolutionary behaviour, see for instance Lax EH . Randall LeVeque E21 and 
Christov and Jordan [23 1. For a detailed exposition on the Godunov method and its exten- 
sions, the reader is referred to [24|. An interesting study on the evolutionary behaviour of 
shocks in elastic nonconductors has been carried out by Fu and Scott [25 , 26 1 using singu- 
lar surface theory together with the method of multiple scales, method of weakly nonlinear 
geometrical optics and the shock fitting method. 

The essential approximation involved in the method of weakly nonlinear geometrical 
optics is that the wave is of small amplitude and is slowly modulated; this means that the 
wave profile changes shape over distances or times much larger than a typical distance or 
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time over which the wave profile itself varies. This method, unlike the method of character- 
istics, is applicable for problems involving several independent variables; indeed, a formal 
and systematic derivation of nonplanar simple wave solutions is possible only in the high 
frequency or geometrical-acoustics limit. Modulated simple wave theories have been de- 
veloped by several authors; papers by Parker l27ll28l and an extensive review by Seymour 
and Mortell [29] are particularly illuminating. For the generalization of the weakly non- 
linear theory in high frequency approximation and for applications of this work, we refer 
the reader to the work carried out by Hunter and Keller |30|, Hunter [31 1, Anile et al [32|, 
Srinivasan and Sharma [33| and Ali and Hunter [34|. 

The other approximation method, namely the simple wave formulation using Riemann 
invariants is useful in the construction of solutions and development of singularities. For 
hyperbolic systems involving source terms or non-planar geometry, the governing equa- 
tions may not admit exact simple wave solutions, but under small wave amplitude assump- 
tion, they may admit asymptotic modulated simple wave solutions. Indeed, hyperbolic 
systems are mathematically tractable using simple waves and Riemann techniques of gas- 
dynamics, see for instance, Zheng [35|, Sterck 11361 and Li et al 071 . Lastly, we consider 
the method of relatively undistorted waves; this applies to solutions of hyperbolic partial 
differential equations, and is based on the same slow modulation approximation used in the 
method of weakly nonlinear geometrical optics. But, unlike the method of weakly nonlin- 
ear geometrical optics, this theory makes no assumption on the magnitude of a disturbance 
(see, Varley and Cumberbatch [38, 39 1 and Seymour and Mortell ll29l ). 

It may be recalled that an analytic description of the evolutionary history of a shock of 
arbitrary strength, both in fluids and solids, still remains an open problem. However, kine- 
matics of a weak shock has been studied in the references cited as above, mostly in solids, 
using singular surface theory and shock fitting methods. The purpose of the present paper 
is to study planar, axially and radially symmetric flows of a polytropic gas behind a weak 
shock, propagating into a uniform state at rest, using four different approximation methods 
and establish their equivalence within their common range of validity. To the best of our 
knowledge, this forms not only the first comprehensive study of one-dimensional planar 
and non-planar weak shocks in ideal gasdynamic flows, but also provides a comparative 
analysis of four seemingly different approaches. The approximation methods employed 
herein for the derivation of asymptotic decay laws for weak shocks and the accompanying 
first order discontinuities are i) the singular surface theory, ii) the method of weakly non- 
linear geometrical optics, iii) the characteristic approach using simple waves as kinematic 
waves, and iv) the method of relatively undistorted waves. 

The paper is organized as follows: after setting up the basic equations and jump con- 
ditions across a shock of arbitrary strength, we devote sections 3 and 4 to the derivation 
of transport equations for the variation of jumps in pressure and its space derivatives; the 
equations are coupled with those involving jumps in higher order space derivatives of pres- 
sure. In section 5, we analyze the pair of coupled transport equations for the shock strength 
and accompanying first and second order discontinuities; the shock strength is assumed to 
be small but the associated first order jump discontinuity may be finite or small. Under 
these assumptions, the transport equations are solved exactly to the leading order once the 
strength of the accompanying second order discontinuity is neglected, and the evolution- 
ary behaviour of both the shock and the associated first order discontinuity is completely 
analysed. Furthermore, the transport equations have also been numerically integrated in 
order to elucidate the results of this section. In section 6, we show how the method of 
weakly nonlinear geometrical optics leads to the same evolution laws for the shocks and 
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first order discontinuities as those obtained in the preceding section. Sections 7 and 8 deal 
with the other two approximation methods, where it is shown that they lead to the same 
evolution laws within their common range of validity. The main results and conclusions 
inferred from this work are given in the last section. 



2. Basic Equations and Jump Conditions across Shocks 

We consider the one dimensional unsteady flow of a polytropic gas behind a planar, 
cylindrical or a spherical wave propagating into a uniform state at rest. The equations gov- 
erning the flow are l40l 

(1) 

Pt + up x + pu x + 3 -pu = 0; u t + uu x + ±p x = 0; p t + up x + pa 2 (u x + ^) = 0, 

where p is the density, u the gas velocity, p the pressure, and a = {"fp/p) 1 ^ 2 the sound 
speed with 7 as the specific heat ratio of the gas; here t stands for time and x the dis- 
tance being either axial in flows with planar (j — 0) geometry, or radial in cylindrically 
symmetric (j = 1) and spherically symmetric (j — 2) flow configurations. 

As system ([T| is quasilinear with each of these equations being a direct consequence of 
the corresponding conservation law, one expects shocks to appear in the flow after a finite 
running length or time. We consider the gas motion containing a shock wave propagating 
into a homogeneous quiescent equilibrium gas with intrinsic velocity U given by U^ 1 = 
s'(x) — -4| where t — s(x) denotes the location of the shock at position x. Across the 
shock front, we have the following Rankine Hugoniot (R-H) conditions, 

(2) U[p] = [pu]; U[pu] = [p + pu 2 }; U[\pu 2 + pe] = [{\pu 2 + pe) u + pu], 

where [/] =/-—/+ denotes the jump in a variable / with / + and /_ being the values 
of / immediately ahead of and immediately behind the shock, and e is the specific internal 
energy given by e = ^ £ . 

In addition, following the singular surface theory, we have the following compatibility 
condition [43], 

(3) ^- = [ft] + U[U 

where -4f denotes the time derivative following the shock front. The medium ahead of the 
shock front is assumed to be uniform and at rest, i.e., p+,p+ and a + are constant, and 
u + = 0. The velocities u, a and U appearing in (JT| and Q are nondimensionalized by 
a + , and the remaining variables p,p, x and t are rendered dimensionless by p + , p+a 2 ^, Xq, 
and respectively; here xq characterizes the reference length of the medium. Thus, the 
variables appearing in ([TJ and (|2| will henceforth be regarded as dimensionless; indeed, in 
the medium ahead of the shock, we have a = a + = 1 and p = p + = 1 . The following 
expressions follow readily from conditions Q, i.e., 

(4) M = m = ™. M = 4r^ 

In view of it follows that, 
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dt (U 2 + l) V dt J V(7 + !)y V dt J 



(5) M = 2U3 (M\ = ( a V . (M\ 



where p = 2 + (7 - l)U 



2 



4 



V D SHARMA AND RAGHAVENDRA VENKATRAMAN 



3. Transport Equation for Shock Strength using Singular Surface 

Theory 

In order to derive the equation for the shock strength, we take jumps in the Euler equa- 
tions ([TJ, and using the conditions ([3]l-(|5]l to obtain 

(6) ^-+A[V x ] + n[u]B = 0, 

where, 

[v] = ( [ P ] M \p] ) tr , [v x ] = ( [ Px ] 

( [«] - U 1 + [p] \ 
A=\ [u]-U i^j , B- 
V 1 + 7 W [u]-U J 

and il = . It may be noticed that the form (|6]l is not only convenient for algebraic 
manipulation, but is also helpful in understanding the pattern of transport equations for 
higher order discontinuities. 

Notice that |6]) contains the unknown quantities [u x ], [p x ] and [pj; eliminating between 
these equations the unknown discontinuities [p x ] and [u x ] by taking a suitable linear com- 
bination, we arrive at the first transport equation governing the shock strength 

(8) ^- = k 11 \ Px ] + k 12 , 

where the coefficients fc n and fc 12 are given by 

2(U 2 - l)p , A(U 2 - 1W / n 

fc n = ^tfttt, — ; — r - ; — > fc i2 = --■ 




U 2 {2p + v) + v 1 U 2 {2p + v) + v V(7 + l 

with v = 2jU 2 + 1 — 7. We note that an immediate consequence of ^ cannot reveal the 
complete history of the evolutionary behaviour of the shocks under consideration because 
of the appearance of an unknown entity [p x ]. It is clear from <(8j that the evolutionary be- 
haviour of shocks at any time t depends not only on the strength of the shock, its curvature 
and the specific heat ratio 7, but also on the pressure gradient immediately behind the 
wave \p x ]. In the following section we determine a transport equation for the unknown 
\p x ] , which represents the effect of disturbances that overtake the shock from behind. 

4. Transport Equation for the Unknown term \p x ] 

As the coupling term [p x ] in ([8| is unknown we need to obtain a transport equation for 
it; to this end, we first note that the [u x ] and [p x ] eliminant of equations |6| along with <(8j 
yields a relation between [p x ] and [p x ]. Similarly the [p x ] and [p x ] eliminant of equations 
(|6]l along with ^ yields a relation between [u x ] and [p x ]. In this section, we determine the 
transport equation for the unknown pressure gradient jump as described in the preceding 
section. We first note that from equations |6| 

where T = [TV/]2x2 is defined by 
(ID 

/ M-(7+l)fclll/ 2 fc 12 (7+l)(l/ 4 -l) 

T_ I vU fen(7((2 M +i/)C/ 2 +i/) 

' I ( 7 +l) 2 t/(A' 2 t/-(7+l)(A'E/ 2 -^)fcll) U(~/+l)((j+l) 2 (U 2 +3)fc 12 +4 M C/([/ 2 -l)n) 

\ vp, 3 2^ 3 
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Under this setting, we now differentiate equations (|TJ, take jumps across the shock, and 
use the shock conditions ^ and Q to obtain the following system of equations 

(12) i^p. +„4[V XX ] +ilV[V x } + tt'[u]B + C = 0, 

where [V] , [V x ] , A, and B have the same meaning as before and 



'XX I 



[V xx ] = ( [pxx] [Uxx] \P: 

( [Vi] (1+[ P }) \ / 2 K][p x ] 

v= o o o , c = -( [Ux p + ^l 

V (1 + 7W) 7(7 -i)M/ V (7 + i)W[pJ 

Eliminating the unknown discontinuity terms [p xx ] and [u xx ] between equations ( fl2| i, and 
using |2]) and (jTOj we arrive at the second transport equation that governs the discontinuity 
\Px], i-e., 

(13) —T^- + k 2l[Pxx] + k 2 2[Px} 2 + k 2 3[Px] + k 2 A = 0, 



d.r 



where 
(14) 



(IP 1) V 

k 2 l = 775 , T) = 



U 2 ' ' 2 M -( 7 + l)C/fc u ' 



~! \ Tl1 V" T FTT); T i,—; \~w)) - f-ttt- 

\ivr[F 22 , (r}vk\\\ 7 + 1 ^ fci 2 dT u . 

' OU-i 12 Z 



/7 4 (7 + l) 2 V 4 M / V fc ii dc/ 



-t/ 2 / 1 J (7+l) 2 U"(7 + l) 



2,(^-1) 

+ ($?) PS + f a^)^ + (7 + 1) (ft?) %T I2 |J 



It may be noticed from ( fT3| l that an analytical description of [p x ] is again obscured by the 
presence of higher order jump discontinuity [p xx ] which is an unknown, whose determina- 
tion requires a repetition of the above procedure; proceeding in this manner we arrive at an 
open system consisting of an infinite set of transport equations for the coupling terms - the 
jumps in higher order derivatives of p. In order to provide a natural closure on the hierar- 
chy of these equations, we consider the coupled system |8]l and ( fTJI and render it tractable 
by making assumptions on the higher order discontinuity [p xx ] (see Fu and Scott[25 1). 



5. Evolution Laws for Weak Shocks 

For a weak shock we may assume that [p] = 0(e), < e << 1, and consider the two 
cases, [p. x ] = 0(1) and [p x ] = 0(e). In the former case, direct perturbation ansatz yields 
a uniformly valid solution to the leading order approximation, while in the second case we 
use the method of multiple scales. We now present these details below. 
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(15) M = [p]+0(6 2 ); [p] = [p]+0(e 2 ); U = ( 1 + [p] ) + 0(e 2 



For a weak shock, it follows from Q that 
and thus the equations ([8) and ( fT3| can be approximated to give 

(16) ^+ j ¥\p\\p*\+m=^ 

%i + ^ww + ^H 2 + iM = o- 

We need to solve these equations assuming that the initial conditions are [p]| x =i = h and 
[Px]x=i = k, where \h\ is assumed to be small and can be taken as e. 

5.1. Case 1: \p x ] = 0(1). When the effect of the first order induced discontinuity or the 
disturbances that overtake the shock from behind is strong, i.e., [p x ] = 0(1), equations 
( fT6l > together with [p] = eY(x), [p x ] = Z(x) and [p xx ] = 0(e q ), q > yield to the leading 
order the following system 

dY 7 + 1 ( j \ _ dZ 7 + 1 ^ 2 



(17) rz HtJ r=0 ' z+ ltJ z ^ ' 

with the initial conditions 

(18) Y(X) = §, Z(l) = k. 

Equation (166) is the well known Bernoulli equation that governs the amplitude [p x ] of 
an acceleration wave; a complete analysis of such an equation modifying and generalizing 
several known results in the literature has been reported in |45| and |46|. We thus con- 
clude that the precursor disturbance that overtakes the shock from behind evolves like an 
acceleration wave at the leading order. System (|T7]>-(p~8|> can be directly integrated to yield 

(19) Y = hl-%ip(x), Z = kr 1 ^{x), 
where ip{x) — x~i and 

(20) 1=1 + ^^fcj(x), J(x) = f ip(s) ds. 



It may be noticed that if k > 0, equations (jT9]» imply the following asymptotic decay laws 
for the shock and the associated first order discontinuity 



1 ( a;" 1 / 2 



x ' , plane 

2 1 1 



(21) ]p]~h[ f— ) ( ji x ~ 3/i > cylindrical, 

' ' J I -7;-i(logx)- 1 /2 ; spherical; 



and 



x , plane , 

(22) [ Px ] - h ( ) { is" 1 , cylindrical, 



7 + 1, 



(xlogx) , spherical, 



as x —> oo; the asymptotic law for the shock is in full agreement with earlier results 
obtained using simple wave theory [40, 44 1. It may be noticed that the decay of the induced 
discontinuity [p. x ] is much faster than the decay of the shock. In order to further elucidate 
the results of this section, the transport equations ( fT6| ) were numerically integrated for the 
planar case for an ideal gas with 7 = 1.4. Figure 1 shows the results of the numerical 
integration. Table 1 further establishes the validity of asymptotic results ( pT| i and ( |22| ) by 
comparing them with results of numerical integration of system ( fT6] ). 
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Numerical Integration of (1E)for h = 0.32, 0.28. gamma = 1.4. 
ID | 1 1 1 1 1 1 1 1 





M : 

W : 
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Figure 1. Numerical solution of system ( fT6| ): Planar Configuration 

h = 0.32, k = 10,7 = 1-4 





h = 0.32, k = 10 


h = 0.32 


,k = 0.28 


Position x 


Error in [p] 


Error in [p. r ] 


Error in [p] 


Error in [pj 


1.476 


4.332e-2 


9.545e-l 


3.374e-2 


6.752e-2 


4.565 


4.827e-3 


4.914e-2 


1.317e-2 


1.885e-2 


7.668 


2.076e-3 


1.593e-2 


7.448e-3 


8.723e-3 


9.563 


I.464e-3 


9.990e-3 


5.675e-3 


6.133e-3 


13.3 


8.752e-4 


5.032e-3 


3.71 le-3 


3.482e-3 


27.95 


2.798e-4 


1.100e-3 


1.373e-3 


9.242e-4 


45.57 


1.332e-4 


4.088e-4 


6.879e-4 


3.678e-4 


65.31 


7.732e-5 


1.979e-4 


4.078e-4 


1.832e-4 


76.04 


6.145e-5 


1.457e-4 


3.271e-4 


1.366e-4 


86.34 


5.075e-5 


1.129e-4 


2.716e-4 


1.065e-4 


96.35 


4.301e-5 


9.054e-5 


2.311e-4 


8.591e-5 


99.95 


4.07e-5 


8.41 le-5 


2.205e-4 


8.072e-5 


100 


4.067e-5 


8.403e-5 


2.190e-4 


7.994e-5 



TABLE 1 . Numerical and asymptotic results for shock strength and ac- 
companying discontinuity 



5.2. Case 2: [p x ] = O(e). When the induced discontinuity [p x ] is weak,i.e., [p x ] = O(e), 
( fl6| ) together with [p] = eY(x), [p x ] = eZ(x) and [p xx ] = 0(e 9 ), where q > 1, yields to 
next order the following system 

<*> f + ^ z + (I) Y - °^ f + ^ <i) z = °^ 

We notice that the straightforward perturbation ansatzes of the forms Y(x) = Y (x) + 
eYi(x) + 0(e 2 ), and Z(x) = Z (x) + eZi(x) + 0(e 2 ) yield secular terms that blow up as 
x — >• oo and render the solutions nonuniformly valid. In order to circumvent this situation, 
we use the method of multiple scales, which requires the introduction of slow variables 
BTIl . Now we take up the following cases separately. 

Planar (j = 0) waves: In this case, we introduce the slow variable, r\ = ex, and introduce 
the formal expansions 

Y = Y (x, r,) + eY 1 (x, V ) + 0(e 2 ), Z = Z (x, V ) + eZ 1 (x, V ) + 0(e 2 ), 
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which on using in (|23) yield 

(24) 



0(1) 
0(e) 



7^ = 0, 
, ox ' 



dZ _ 



0. 



dx + 9»7 + 2 Z - u - 

with the following initial conditions 

(25) y (l) = Zo(l) = f, 
where \h\ and |fc| are small. The first pair of equations in ( p4| ) yield 

(26) F = /(r?), Z =g(n). 

Furthermore, in order to avoid secular terms from the equations at this order 0(e), we have 

(27) f( v ) + ^.fivMv) = o, g'(v) + ~g 2 (v) = o. 

Finally solving d27| for k > in conjunction with the initial conditions (p5|, we obtain 



(28) [p] ~ ft 



1/2 



[Pa 



as x — )• oo. 



.(7 + W " ' (7 + 1)' 

In order to further elucidate the results of this section, the transport equations ( p"6] ) were 



Numerical Integration of (16) with 11=0.32, k = D.28 gamma = 1.4 




100 200 300 400 500 BOO 700 BOO 900 1000 



FIGURE 2. Numerical solution of system ( fTo*] ): Planar Configuration 
h = 0.32, fc = 0.28,7 = 1.4 



numerically integrated for an ideal gas with 7 = 1.4 Figure 2 shows the results of the 
numerical integration of ( |16) , Table 1 further establishes the validity of the asymptotic 
results in |28) with the results of numerical integration of ( fl~6*| ). It may also be noticed that 
in this case, the shock decays faster as compared to the previous case when k = 0(1); this 
has been shown in Figure 3. Finally, it is interesting from ( |28| i and |22) that the asymptotic 
behaviour of the precursor disturbance [p x ] is independent of k, this has been illustrated in 
Figure 4 where the two curves are indistinguishable in the asymptotic limit. 
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Figure 3 . Comparison of shock strength decay in case 1 and case 2 above 





1 --, 1, t= III 

[p y ] with k = 0.28 . 
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Figure 4. Comparison of [p x ] in case 1 and case 2 above 



Cylindrical (j = 1) wave: In this case, the system of equations ( |23] l with the slow variable 
f] = e^/x yields the following hierarchy of partial differential equations 



(29) 



O(e): ^ + ^f l + ^^o + (i>i = 0, 

dx drj 



^o 2 + (^)^i=0, 



together with the initial conditions ((25). We find that to the leading order Yq = A(rj)x 1 ' 2 
and Zo = B^x^ 1 ^ 2 , where A and B satisfy the differential equations 



A 1 + 1 -^^AB = 0, B' + (7 + l)B 2 = 0. 
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Here prime denotes differentiation with respect to the variable 77; solving these equations 
and using ( p5| , it follows that for k > 0, the shock and induced discontinuity decay asymp- 
totically as 

(30) M~ /t ( fc ( 7+ i) j x " 3/4 ' W~^i) aj_1 asx^oo. 

Spherical (j = 2) waves: We repeat the above procedure with the slow variable rj = 
e log x, appropriate to this case, and find that for k > 0, the decay behaviour to the leading 
order is 
(31) 

/ 2 \ 2 
[p] ~ ft. [ k , +1 \ J x-^logxy 1 / 2 , \p x ] ~ aT^logz)' 1 as x -> 00. 

6. Weakly Nonlinear Geometrical Optics 

In this section, we use the theory of weakly nonlinear geometrical optics, which allows 
waves to propagate on different families of characteristics, to treat shock formation, and 
to include interaction of the flow with the shock. This is an asymptotic method whose 
objective is to understand the laws governing the propagation and interaction of small 
amplitude high frequency waves QUI . Here we show that this method leads to the same 
evolution laws for shocks and first order discontinuities as those obtained in the preceding 
section. 

We assume that the disturbance is of small amplitude as in the previous section; the 
parameter e used in the previous section suggests to form an asymptotic analysis. We take 
(x, t, 9) as independent variables where 9 — is the fast variable with cf>(x, t) as the 

phase function to be determined, and = constant is the characteristic wavelet which left 
the boundary x = 1 at 6 = -. We now seek small amplitude high frequency solutions to 
the system ([T} of the form 

(32a) p=l + ep (1) {x, t, 9) + e 2 'p™ (x, t, 9) + 0(e 3 ), 

(32b) u = eu (1) [x, t, 9) + eV 2) (x, t, 9) + 0(e 3 ), 

(32c) p = i + ep {1) {x,t,9) + e 2 p (2 \x,t,9) + 0(e 3 ). 

where p^ 1 ', vP-\p^ etc are regular and bounded functions of their arguments together 
with their derivatives as e — > 0; using ([32} in ([1} along with 

?l\ = dfi 1 , 91. df\ = 0/i ,l 6 ?l 

dx \t dx^ e e 9x d0' dV* dt^' e e 9t dQ ' 

and equating to zero the coefficients of e 4 , i = 0, 1, we arrive at the following hierarchy of 
partial differential equations: 

(33a) 0(1) : <f> x u^ + cj> t p^ = 0; 

(frxP^ + <th^e = 0; 

<j>xV>g + <PtP { e ] = 0, 
(33b) O(e): Pt (1) + <i> t pf + <p x ({p (1) u {1) )e + u™) + u<?> + ^ = 0; 

+ <t> t uf + + pf " P {1) vf ] ) +Pl l) = 0; 
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Equations q33p admit nontrivial solutions if <p satisfies the eikonal equation 



which shows that the phase function <fi(x,t) = constant represents a propagating wave front 
of the leading order system associated with ([TJ. Indeed, each of the phase speeds, — j 1 , 
which can take the values —1,0 and 1, gives the speed of propagation of <fi along normal 
to the wave front. It follows from ( |33a[ ) that along the forward facing wave, — = 1, we 
have 

(34) =uW =pW =a(x,t,6), 

where the wave amplitude a remains undetermined to this order of approximation. It may 
be remarked here that the propagation along the other two families of characteristics can 
be handled in a similar manner. Using ([34} in the p^ 1 , and pjp eliminant of ([33b} we 
get the following transport equation for the wave amplitude: 

m\ , a f^ + 1 \ . 3 a n 

(35) a t + a x - <p t I - aag -" 



2 J 2x 

Furthermore, if we impose the condition that <f>\ x =i = t, then — ^- = 1 yields 
(36) <j) = t-{x-l). 

Accordingly, we may solve ([35} to obtain 



I 

(37) ax 2 = constant , 
along the family of curves given by 

ns\ dx 1 d6 7 + 1 

(38) -dt =h itt=-^r a - 

Let these curves (characteristic wavelets/or rays) be marked by the characteristic variable 
r, which is chosen such that t = r at the boundary x = 1; then integrating ( [38} and using 
the definition of 9, we find that the equation of any wavelet, r = constant , is 

(39) t-r = x- 1- ^4r~ u ^ J(x), 



where 




plane 



(40) J(x)= I r~2dr={ 2{y/x-l), cylindrical; 

spherical. 

In writing ( [39} , we have replaced eu^ by u, as it is true in the leading order approximation. 
Also, equation ( [37} can be written as 

i 

(41) UX2 = v (t), 

where v(t) is the value of u on the boundary x = 1. Equations ( [39} and ( |4l} constitute 
the desired modulated simple wave solution; infact, it follows that along a positive charac- 
teristic that intersects the boundary at time r, we have u = v(t), and the equation of the 
characteristic is 

(42) t - t = x - 1 - ^yiw(r) JO). 
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Here, we envisage a weak shock wave that propagates into an undisturbed region with 
u = ahead of the shock. Thus, it follows that across the shock 

(43) [u] = x~iv(r-), 

where r_ is the wavelet immediately behind the shock. Also, if the shock path is t = s(x), 
then the shock conditions Q implies that to the leading order approximation 

(44) — = s' (x) = U- 1 ~ 1 - L — v(t-)x~2, 
ax 4 

and the characteristic equation ( [39] > gives 

(45) s(z)-t_ =x- 1- 2^u(t_)J(x). 

Equations ( |44| and ( |4"5j ) determine the shock path £ = s(x) by elimintating the parame- 
ter r_. Further relations for the shock speed U follow from ( |45) ; this together with 
implies that on the shock 

7 + 1/ x -i 7+1./ nt/ x 

(46) V U(T - )S 2 + — V ^ X) ^ = 
which yields on integration 

(47) ^±l v i( T _)j( x ) = f v ( s )ds. 

4 Jo 

Also, since the characteristics r + = ahead of the shock and the characteristic t_ behind 
the shock meet at the shock at the same point and at the same time and since v(t+) = 0, it 
follows from d42l» that 



'>o\v(t-)J v(0) 2 

It may be noticed that ( |48| i also follows from (j47j in the limit t_ — > 0, showing thereby 
that the integral in ( |47| i is bounded at large distances (i.e., x — > oo.) Thus, if the initial data 
is such that v(t) — > as r — > To, then some of the characteristics in the neighborhood of 
the zero wavelet do not overtake the shock from behind, and the asymptotic behaviour as 

x — > oo is 

/ 4& N 1/2 

( 49 ) u(t_) ~ 7 — T . 

where | JJ" v(r) dr\ <b. Relation ( |49| ) together with ( |47| ) implies that 

_j ( 46 \ 1/2 _i 

(50) M= U -= V ^ X -~ 2 ~{ {7+1)J{x) ) x ~~ 2 ' 

In other words, the shock decays like 

x" 1 / 2 , planar; 

(51) [u] ~ { x-*^, ^ cylindrical ; 

" 1 (log x) ~ 2 , spherical , 

as x — >• oo; this result is in agreement with that obtained in the preceding section. Further- 
more, it follows from |42) that between the shock and the limiting characteristic r = To, 
we have 

(52) x ~ t - r + 1 ^v(t) J(x), 
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so that 

a; -1 , plane ; 

(53) u = v(t)x~z ~ — T~\( x — ( r — r o)) ■( s^ -1 ) cylinder; 

(xlogx)" 1 , spherical. 



^7 + 1' 



The asymptotic behaviour of the first order discontinuity at the shock then follows from 
( |5"3j ) by differentiating it partially with respect to x, and then evaluating at the shock, 
namely 
(54) 



2 



7' 



x 1 , plane ; 

-{K{x)+{x-s{x)+t q )K'{x)}\ K{x) = I ix" 1 , cylinder; 

(xlogx) -1 , spherical. 

Evaluating ( [52] ) on the shock and then using it in ( |56| l above, we find that to the leading 
order approximation 

x _1 , plane; 
(55) [u x ] ~ — — { 5 s-1 j cylindrical; 



2 

' I (xlogx) -1 , spherical, 



as i — > co, which agrees with the decay behaviour outlined in the preceding section; 
indeed, it can be easily verified that in the weak shock limit [u x ] ~ [p x ] to the leading 
order. 

7. Simple Wave Formalism using Riemann Invariants 
System ([TJ can be written using vector matrix notation as 
(56) u t + Au x + B = 0, 

where u = ( p u p)* r ,_B=( jpa2 " )' r and A is a 3 x 3 matrix that can 

be read off by an inspection of ([TJ. Since the eigenvalues of A are u ± a and «, and the 
corresponding left eigenvectors are ( ±ap 1 ) and (a 2 — 1 J , respectively, it 
follows that equations d56b can be written in the characteristic form as 



,™ dp , du jpua 2 dx 

(57) -£ ± pa— + ^ = along C± : — = u ± a, 

OP 2 dp n , n dx 

(58) — -a 2 -£=0 along P : — = u, 

or at at 

where C± represent characteristics moving with velocity ±a relative to the velocity u, and 

P the particle path. Since p = p(p, S) is a given function, where S denotes the entropy, 

equation |58| shows that the entropy remains constant along each particle path 4| = u. 

Further, as the jump in entropy at a shock is of the third order of smallness relative to the 

discontinuity in pressure, equation |58| implies that for a weak or even moderate strength 

shock, S = So (constant) is a good approximation everywhere, and consequently equations 

(|57j can be written as 

(59) — — + -au = along : — = u ± a, 

at x at 

where = i u ) am tr| e Riemann invariants, and the differential operators on 

are just the operators of differentiation along the characteristics C ± . The flow described 
by the solution of ([59j> is a simple wave; in the plane (j = 0) case, equation (57b) is 
dispensed with quickly by arguing that r~ is constant everywhere. For a nonplanar case, 
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the situation is different. As the jump in r~ across a shock is also of the third order of 
smallness relative to the jump [p], and the integral of the last term in (57b) taken along 
C~ is small keeping in view that the simple wave argument applies in general near the 
wave front of any disturbance propagating into a uniform state at rest, it follows that the 
constancy of entropy S and the Riemann invariant r~, i.e., S = Sq and r~ = -^j are a 
good approximation to the first order, and thus the first equation in ( |57| ) provides a single 
first order equation for u, namely, 

(60) ut + (l + "^uj u x + J -u(l + ^y^t) = 0. 



Equation (|60j can be solved exactly to yield 

2 

/ 7-I \ 7-1 £ 

(61) ii 1 + ^— ii =v(t)x 2. 



where r is the characteristic variable to be determined from ^ = (l + -2^) 



7+1 „ 
2 ' 

Thus, a uniformly valid approximation at the leading order is 



_l 

(62) u = v{t)x 2, 

t = T+(x-l)-^v( T )J(x), 

which are exactly the same as (39)-(41), and hence the evolutionary behaviour of shock and 
the first order discontinuity turn out to be the same as discussed in the preceding sections. 

8. Method of Relatively Undistorted Waves 



This method proposed by Varley and Cumberbatch [ 38 . 39 1 is based on a scheme of 
successive approximations to a system of hyperbolic equations. The theory does not make 
any assumptions on the amplitude of the disturbance and yields exact results for simple 
waves, acceleration waves, and thus for the formation of shock fronts. Here, we show that 
this method yields the same results as those derived in section 7 for the evolution of shock 
strength and the accompanying first order discontinuity. 

Under the transformation of coordinates (x,t) (x,t), where r(x,t) is a wavelet 
defined by t = t(x,r), and any vector valued function u(x,t) transforms as u(x,t) H > 
U(x, r), system |56) may be written as 

(63) (A-t x 1 l)u t =t- 1 {AU x + B), 

where 1 is the identity matrix of order 3. A solution vector u is said to be relatively undis- 
torted if 

(64) ||U X || « K||, 

implying thereby that u x ~ t x u t . Equations ( |63j ) and (|64j> are compatible if ||B|| = 
0(1) ll^Ua, || , while t^ 1 is an eigenvalue of A: this implies that the wavelets are the char- 
acteristic curves of ( |56l >. Subsequently it follows from ( |63] l that U must satisfy the com- 
patibility condition 

(65) HAU X + B) = 0, 

where L = ( pa 1 ) is the left eigenvector of A corresponding to the eigenvalue 
t^ 1 = u + a. We look for the solution of ( |56| ) in the region x > 1, where the motion 
associated with the eigen-mode t x x = u + a is excited by the boundary condition 

(66) u(1,t)=v(t), 
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where v is a smooth bounded function, i.e., |y| = 0(1). We let r = constant denote the 
characteristic curve defined by t^ 1 = u + a, and parametrize r by defining r = t at x = 1. 
Thus, equation ( |o*3j ) to a first approximation yields 

(67) (A-t- 1 2)XJ T = 1 

which together with ( |65) and t~ x = u + a enables us to determine the unknowns U and 
t(x, t) . Equation (|67) implies that U is collinear to the right eigenvector R = ( p a pa 2 
of A corresponding to the eigenvalue u + a. For a wave moving into an undisturbed region, 
equation |67) integrates to give 

2 2 7 
/ 7-I \ 7-1 If 7-1 \ 7-1 7-I 

(68) p = ( 1 + -L-^-u J , p = - f 1 + l^-tt J , a = 1 + 



Equations (68) hold throughout the region at any a: on any wavelet r = constant; these 
may be used in ( |6"5) to yield the following transport equation for u, i.e., 

1 _l T + 1 ^ j /^i j. T _ M n 
1 H u w T H u 1 H u = U, 

2 / 2x I 2 / 



which on integration yields 



2 

7 - 1 \ 7-1 ri 



2 



u(r)a; 2 ; 



where r is the characteristic variable to be determined from t x 1 = (1 + " ~ 

(l — ~2 L u) ; this indeed, is precisely the solution (61). Thus, the relatively undistorted wave 
approximation is a consequence of the simple wave formalism using Riemann invariants. 

At this juncture, it is worthwhile to remark that the Chester-Chisnell-Whitham (CCW) 
approximation, which consists in applying the shock conditions Q to the exact character- 
istic equation ( |57] > along C+, yields the following transport equation for the shock 

(69) l7^T^ + x =0 ' 

where -4- denotes the derivative along C+ : = (u + a)| s hock' an< ^ = (l + 

2(^) 1 ^ 2 + J7 _2 ){1 + ^xrj2 }■ It is indeed interesting to note that the transport equation 
(8) for the shock, in view of relation (ffl, assumes the following form by neglecting the 
precursor disturbances from the rearward flow (i.e., on setting [p x ] = 0) 

where G(U) = (7 + 1){2U 2 ^- 1 + p^^U 2 + 1)}; it may be noticed that (|69j and (70), 
except for the functional forms of g(U) and G(U) and the differential operators, share a 
close structural resemblance. Indeed, it is remarkable that for weak shocks, U ~ 1, both 
G(J7) and g(?7) ~ 4, and the evolutionary behaviour governed by ( |69) and (70) coincide. 

9. Results and Conclusion 

Except some impirical methods, which have been developed in the past, no analytical 
method exists for describing the evolutionary behaviour of a shock wave without limiting 
its strength. In the present paper, we generate a completely intrinsic description of plane, 
cylindrical and spherical shock waves of weak or even moderate strength, propagating into 
a uniform flow of a polytropic gas at rest. We use four different approximation methods 
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to study the asymptotic evolutionary behaviour of shocks and the associated first order 
jump discontinuity, which may be interpreted as the disturbance that catches up with the 
shock from behind. Singular surface theory yields a set of coupled equations for the shock 
strength and the higher order jump discontinuities which couple the shock motion with the 
rearward flow. It is noticed that the transport equations for the shock strength [p] and the 
first order jump discontinuity [p x ] are nonlinear in [p] and [p x ] respectively; however, the 
transport equations for [p xx ] is linear in [p xx ], and this is true for all other higher order 
jump discontinuities. The transport equation for the shock strength shows that the evolu- 
tionary behaviour of the shock at any instant is influenced by the shock curvature and the 
first order jump discontinuity, representing the effect of precursor disturbances that over- 
take the shock from behind. This transport equation, in some sense, generalizes the CCW 
approximation [40], which although does not account for the effects of the disturbances 
that overtake the shock from behind, yet its success is unexpectedly remarkable for a par- 
ticular class of problems. We show in the preceding section that the transport equation for 
the shock strength, in the limit of a weak shock and vanishing [p x ], reduces exactly to the 
one obtained by the CCW approximation. The truncated system of transport equations for 
[p] and [p x ] is solved for uniformly valid solutions using a singular perturbation technique 
fiTTl that eliminates secular terms in the perturbed solution; the shock strength is assumed 
small, but the accompanying first order jump discontinuity may be finite (i.e., 0(1)) or 
small. It is shown that the precursor disturbances evolve to the leading order like an ac- 
celeration wave. Indeed, the decay of the induced discontinuity is much faster than the 
decay of the shock; it is noticed that for the case [p x ] = 0(e), the shock decays faster as 
compared to the case when [p x ] = 0(1). Asymptotic results for the decay of weak shocks 
fully agree with those given in fiOl |44l . It is shown that the approximation methods of 
weakly nonlinear geometrical optics, simple waves involving Riemann invariants l42l . and 
the relatively undistorted wave approximation yield the same results within their common 
range of validity. 

In conclusion, we observe that an analytical description of the evolutionary behaviour 
of shocks of arbitrary strength in an ideal gas remains an open problem because the trun- 
cation approximation, used in rendering an infinite hierarchy of transport equations as a 
closed finite system, may not yield a uniformly valid solution for all running lengths. The 
present work envisages four different methods to study the evolutionary behaviour of weak 
shocks and the accompanying first order discontinuity; a thorough comparative analysis 
of these methods has been presented, which to the best our knowledge has not been done 
earlier. Indeed, the singular surface theory gives much more powerful results then does the 
method of characteristics that forms the basis of rest of the approximation methods used 
herein. The truncation approximation leading to the closed system consisting of a coupled 
pair of transport equations for [p] and [p x ] can be regarded as a good approximation for the 
hierarchy of the system governing shock evolution. Further, as the relatively undistorted 
wave approximation, and hence the method of weakly nonlinear geometrical optics, which 
satisfies the condition needed for an approximation to be relatively undistorted, lead to the 
first order the same evolution laws for the shock and the associated first order discontinu- 
ity, it transpires that both these approximation methods are a consequence of the simple 
wave formalism using Riemann invariants. However, the relatively undistorted approxi- 
mation makes no assumption on the amplitude of the disturbance, unlike the later which 
is a manifestation of the method of multiple scales. Lastly, a comparison with the CCW 
approximation reveals that the singular surface theory yields the shock-transport equation 
which can be regarded as a generalization of the CCW approximation. 
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